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functions of one single variable x = t/^o- These predictions are confirmed by accurate 
numerical data, both in the normal and in the anomalous case. 
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1 INTRODUCTION 



The Anderson localisation in a random potential is now well understood, at least 
as far as static or spectral properties are concerned [1]. In the one-dimensional case, 
all eigenstates are exponentially localised, with the localisation length C{E) depending on 
energy E [2, 3]. Dynamical aspects, concerning mostly the spreading of a wave packet, have 
also recently attracted considerable interest [4-6], especially in connection with diffusion 
in driven quantum systems and random band matrices [7-9], and more recently with the 
problem of two interacting particles [10]. 

In this paper, we emphasise a novel and striking feature of the dynamics of the one- 
dimensional tight-binding Anderson model. The q-th moment of the local probability 
density, or dynamical participation number of order q, to be defined in eq. (1.4), exhibits 
anomalous scaling and bifractal behaviour, with a non-trivial exponent T{q) for q > 2. 
This phenomenon seems to have been entirely overlooked so far. It will be shown to take 
place in the absence of disorder, and in the ballistic regime, in the localised regime, and 
throughout the crossover between them. 

To be more specific, we investigate the time-dependent wavefunction of a tight-binding 
electron in one dimension, which obeys 

i^^^ = i^n+l{t) + ll^n-lit) + Vn^t), (1.1) 

in reduced units, such that lengths are measured in units of the lattice spacing, and energies 
and inverse times in units of the hopping integral. The diagonal site potentials {fn} are 
independent random variables, drawn from a common distribution. We choose the initial 
condition of a particle sitting at the site n = at the origin of times: 

V'n(0) = 5n,0. (1.2) 

We characterise the spreading of the wave packet by the following quantities. 
• Moments of the position of the particle: 



M,(t) = (|nr) = 5^|n|^P„(0. (1.3) 

n 

• Moments of the probability density (dynamical participation numbers): 



s,{t) = Y.{Pn{t)y. (1.4) 

n 

In the above definitions, the probability density reads 

Pn{t) = \^n{t)\\ (1.5) 

and the index q is any real positive number, not necessarily an integer. The bar denotes an 
average over the disorder, i.e., over the distribution of the random site potentials {vn}- The 
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conservation of the norm of the wavefunction ensures that Mo{t) = Si{t) = 1 at all times 
t > 0. The most commonly considered quantities in the literature are the mean squared 
position M2{t) and the participation number (or inverse participation ratio) 5*2 (t) [11-13]. 

Scaling properties are known to hold in the weak-disorder regime, namely for a small 
enough random potential. Assuming that the site potentials have zero average, it is suffi- 
cient to characterise their distribution by its variance: 



To lowest order in perturbation theory, the localisation length is maximal in the vicinity 
of the band centre, where it scales as [1-3] 



This is the characteristic length scale where the localisation phenomenon takes place. The 
absence of an intermediate diffusive regime between the ballistic one {t -C ^o) and the 
localised one {t » ^o) is a peculiarity of the one-dimensional situation, where there is no 
fundamental difference between the mean free path and the localisation length [14-18]. A 
more detailed description of the localisation length, including its anomalous scaling near 
band edges, will be recalled in section 3.1. 

We shall be interested in the long-time behaviour of the moments Mq{t) and Sq{t) 
defined above, in the weak-disorder regime. Scaling properties can be expected in this 
situation, where both characteristic length scales t and are simultaneously large. The 
setup of this paper is as follows. Section 2 contains a detailed analytical investigation of 
the problem in the absence of disorder. The wavefunction is given by a Bessel function, 
which exhibits three different kinds of asymptotic behaviour in the (n, t) plane. As a 
consequence, the moments Sq{t) of the probability density exhibit anomalous growth with 
a scaling exponent T{q) for q > 2. This behaviour, interpreted as bifractality, is expected 
to hold in the ballistic regime, namely for t -C ^o, where disorder will have hardly any effect 
on the ballistic motion of the particle in the absence of disorder. In section 3 we show on 
a heuristic basis that the same kind of anomalous scaling, with the same exponent T{q), 
also takes place in the localised regime {t 3> ^o): but with the asymptotic width of the 
wave packet, of order replacing time. Finally, we demonstrate that all the quantities 
under consideration exhibit scaling behaviour throughout the crossover between the free 
(ballistic) and the localised (insulating) regimes, involving universal scaling functions of 
the variable x = t/^Q. This prediction is confirmed by accurate numerical data. Section 4 
contains a brief discussion. 

2 ANALYTICAL RESULTS IN THE ABSENCE OF DISORDER 

2.1 Description of the wavefunction 

In the absence of disorder, the dynamics of the Anderson model can be investigated 
analytically. Let us denote the various quantities with the superscript (0) in this limiting 
case. The stationary tight-binding equation reads 



t,2 = < 1. 



(1.6) 




(1.7) 




n-l - 



(0) _ 



(2.1) 



3 



Its eigenfunctions are the plane waves ipn^ = e'"^, where momentum p is related to energy 
E by the dispersion relation 

E = 2cosp. (2.2) 
We thus obtain an explicit expression for the time-dependent wavefunction, 

i(;l^\t) = j ^e^^^'-2'*^°^P = i-^Jn(2t), (2.3) 

where the momentum integral runs over the Brillouin zone B = [—it, tt] . The probability 
density at site n, 

Pi'\t) = \i;^^\t)\' = {U2t))', (2.4) 

is thus the square of the Bessel function Jn{2t), whose argument is proportional to time, 
while its order is the number of the site, i.e., the distance travelled by the particle from its 
starting point. 

It turns out that three regions have to be considered in the (n, t) plane, where the 
Bessel function admits different kinds of asymptotic behaviour. The existence of these 
regions can be explained by the following semi-classical argument. The dispersion rela- 
tion (2.2) corresponds to the group velocity v — dE/dp = — 2sinp. As a consequence, for 
a wave packet initially peaked in momentum space around some mean p^, the centre of 
mass will move according to the semi-classical law 

(n) pa — 2tsinp*. (2.5) 

We thus expect an allowed region (|n| < 2t), separated from a forbidden region (|n| > 2t) 
by sharp fronts located at n f=:i ±2t. 

The above heuristic picture can be made quantitative by means of the following asymp- 
totic formulas in the theory of Bessel functions [19]. As a matter of fact, the derivation 
of these formulas relies on the method of steepest descent, with the saddle-point equation 
coinciding with eq. (2.5). 

• Allowed region {\n\ < 2t): 

For n > 0, and with the notation 

n = 2tsmp (0<p<7r/2), (2.6) 

in agreement with the semi- classical law (2.5), we have 

p(0) ^ sin^ (n(p + cotp-7r/2) + 7r/4) 

^ ^ ^ TTtCOSp ' ^ ' ^ 

The probability density is thus the product of 1/t by a rapidly oscillating amplitude, as 
expected in an allowed region. 
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• Forbidden region (|n| > 2t): 

For n > 0, and with the notation 

n = 2tcosh^ (6'>0), (2.8) 

we have 

.r,^ , . exp ( — 2n(9 — tanh6')) 

The probabihty density thus decays exponentially, as expected in a forbidden region. 

• Transition region (|n| ~ 2t): 

The transition region, corresponding to the ballistic fronts, turns out to extend over 
a spatial range of order t^/^. For n > 0, and with the notation 

n^2t + t^/^z, (2.10) 

the probability density is approximated as 

p(«)(t)^t-2/3(Ai(^))', (2.11) 

with Ai(^) being the Airy function. The asymptotic behaviour of this function, namely 

Ai(.)J-"^l^l"''^'''(iW^'^ + j) (2.12) 
\|7r-i/2^-V4 exp (-1^3/2) (^^+oo), 

respectively matches eqs. (2.7) and (2.9). 

Throughout the transition region, the probability density is the product of by an 
amplitude which oscillates toward the allowed region {z — > —oo), and falls off exponentially 
toward the forbidden region {z — > +oo). 

The existence of these three regions is illustrated in Figure 1, showing the probability 
density Pn\t) against site number n, for a time t = 100. 

2.2 Moments of the position 

We now turn to the analysis of the moments Mg*^'* (t) of the position of the particle, in 
the absence of disorder. These quantities are quadratic forms in Bessel functions, so that 
their analysis is rather easy, at least for even integer values of the index: q = 2k. We have 
indeed 



^ f ^ f V_g2it(cosp'-cosp) ^^2fegin(p-p') 



2n In 27r 



n 



(2.13) 



L 



_-2itcosp ( \ 2itcosp 



2tt^ V V 



The integrand of the last expression can be expanded as a trigonometric polynomial. As 
a result, the moments M^^^t) are even polynomials of t, with positive integer coefficients. 
We have MQ^\t) = 1, as expected, and 

M^^\t) = 2^2, M^^\t) = + 2^2, M^^\t) = 20t^ + 30t^ + 2t'^, etc. (2.14) 

The long-time behaviour of the moments is obtained by letting all the derivatives act 
on the exponential in the last expression of eq. (2.13). We thus obtain 

^2?W~02fct^'' (t>l), (2.15) 

with 

a.=/J(4sin^p)' = p. (2.16) 
This estimate can be shown to hold true for any real g > 0, namely 

M^°)(t)«aqt« (t>l), (2.17) 

with 

2<i r (^) 

«^ = ^ (2.18) 
where r{z) denotes Euler's gamma function. 

2.3 Moments of the probability density (participation numbers) 

The analysis of the long-time behaviour of the moments Sq'\t) of the probability 
density is slightly more involved. Indeed these quantities are highly non-linear fimctionals 
of the wavefunction. As the index q gets larger, they are more and more sensitive to large 
values of the wavefunction, whereas the moments Mq'^\t) of the position are not. 

The asymptotic expressions (2.7), (2.9), (2.11) show that the probability density 
pji^\t) scales as 1/t in the allowed region (bulk of the wavefunction), over an extent of 
order t sites, while it scales as in the transition region (fronts of the wavefunction), 

over an extent of order t^^^ sites, and it is negligible in the forbidden region (tails of the 
wavefunction). Hence the bulk has a normal contribution to the moment Sq^\t), scaling as 
t~ii~'^\ while the anomalous contribution of the fronts scales as This analysis 

therefore predicts the power-law behaviour 

S^^'Xt) ^ bgt-^^'i^ (t>l), (2.19) 

with r(g) being the smaller of both exponents, namely 

q — 1 for q < 2 (normal), 

r{q) = { 2q-l (2.20) 
— - — for > 2 (anomalous). 
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This prediction is summarised in Table 1. It can be made more quantitative by the 
following analysis, yielding the value of the prefactor bq of the formula (2.19) in either case. 

• Normal regime (g < 2): 

In the normal regime, the moment Sq^^ (t) is dominated by the bulk of the wavefunc- 
tion, corresponding to the allowed region. The prefactor of eq. (2.19) can be estimated 
from the expression (2.7), assuming that the arguments of the sine functions are uniformly 
distributed, and transforming the sum over n into an integral over p. We thus obtain after 
some algebra Sq^\t) ^ hqt~^^~^\ in agreement with eq. (2.20), with 

2r(V) rfa + i) 
- ^« r (V) ni + 1) ' ' 

• Anomalous regime {q > 2): 

In the anomalous regime, the moment 5'^'^'' (t) is dominated by the fronts of the wave- 
function, corresponding to the transition region. The prefactor of this moment can then be 
estimated by using the expression (2.11), and transforming the sum over n into an integral 
over z. The outcome again agrees with eq. (2.20), and yields 

bq = 2 dz\Ai{z)f'^ (g>2). (2.22) 



We have in particular 

h = 0.073214. (2.23) 

• Marginal case {q = 2): 

This borderline case corresponds to the usual participation number S2{t). The expo- 
nents of the contributions of the bulk and of the fronts have the common value t(2) = 1. 
It is worth noticing that the prefactor bq diverges as bq ~ 3/(277^10" — 2|) as g — > 2 from 
both sides. This can be checked for q < 2 directly from eq. (2.21), and for q > 2 from the 
behaviour (2.12) of the Airy function as 2 — > — cxo. 

The moment S^2^\t) turns out to exhibit a logarithmic correction to its leading 1/t 
behaviour, which can be analysed by means of the Mellin transformation. The function 
S2^\t) and its Mellin transform m{s) are related by 

poo I* 1 

m{s) = t^-'si'\t)dt, Si'\t) = J ^^t-'m{s), (2.24) 

for Re s positive and small enough. We have 

si'\t)^^{U2t))' 



/■dpi r 

Jb 27r "Jb 



g2it(cOSpi-|-COSp2— COSps— COSP4) giw(p3+P4— Pi — P2) 



27r 

n 



f du f dv f 
Jb2^Jb^ Jb 



(2.25) 

V ' 

27r i5(pi -I-P2 -P3 -P4) 



d^ p— Sit sin u sin v sin w 

s27r' 
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where the last expression has been obtained by the change of variables 

w = (Pi -P2 +P3 -P4)/4, V = {pi-p2-ps+P4)/4:, ^2 26) 

w ^ {n-pi- p2)/2 = (tt - p3 - p4)/2. 

The product structure of the last expression of eq. (2.25) makes it suitable to the closed- 
form evaluation of the Mellin transform m{s). We indeed obtain, for < Res < 1, 



/ \ -r^/ \ [ du [ dv [ dw ... . . . 

m{s) = T[s) 1 — / — / — (Sismttsmvsmtu) 
Jb 27r Jb 27r Jb 27r 



The long-time behaviour of S^\t) is given by the double-pole singularity of m{s) at s = 1, 
of the form 

where 7^ denotes Euler's constant, hence 

We have thus derived both the prefactor and the finite part of the logarithmic correction 
of the dynamical participation number 82^^ (t) in the absence of disorder. The finite part 
is a surprisingly large number: writing the numerator of eq. (2.29) as ln(t/to), we have 
1/to = 64exp(7s) = 113.989. 

2.4 Interpretation: bifractality of the probability density 

The scaling law (2.19), (2.20) for the moments Sq^\t) of the probability density in 
the absence of disorder, with its two branches of exponent T{q), can be interpreted within 
the multifractal formalism [20]. Indeed, the wavefunction takes appreciable values over a 
number of lattice sites of order t. As a consequence, 1/t can be viewed as a short-distance 
cutoff in the definition (1.4) of the moments Sq{t). 

The scaling exponent T(g) can be interpreted in terms of generalised (Renyi) dimen- 
sions Dq of the local probability density. The relation T(g) = {q — l)Dq yields 

for g < 2 (normal), 

(2.30) 

for g > 2 (anomalous). 

The probability density can be alternatively characterised by a multifractal spectrum f{a), 
which is the Legendre transform of the exponent T(g), according to 




The expression (2.20) yields the following results. The normal branch T{q) = q — 1, i.e., 
Dq = 1, for q < 2 yields the point (a = 1, / = 1), corresponding to the normal scaling 
of the bulk of the wavefunction, while the anomalous branch T{q) = {2q — l)/3 for q > 2 
yields the point {a — 2/3, f — 1/3), corresponding to the anomalous scaling of the fronts 
of the wavefunction. 

These results are summarised in Table 1. We propose to call bifractality such a scaling 
behaviour, with a normal and an anomalous component. 

3 SCALING ANALYSIS IN THE GENERAL CASE 

3.1 A reminder on band-edge anomalous scaling 

We first recall some results on the scaling behaviour of the localisation length in the 
presence of a weak diagonal disorder. Inside the band of the pure system, characterised by 
the dispersion relation (2.2), and in the weak-disorder regime (d^ <^ 1), the localisation 
length scales as [1-3] 

8 sin^ P ^ 9 /X 
C^—^^^Qsm^p. (3.1) 

This leading-order perturbative prediction vanishes as p — > or p — > tt, corresponding 
to the band edges, namely E — > ±2. This observation suggests that the localisation phe- 
nomenon has something special near band edges. This effect has been initially investigated 
by Derrida and Gardner [21], who indeed demonstrated the presence of anomalous scaling 
in the localisation length ^{E) and the density of states p{E). These quantities behave 
near the upper band edge {E —>■ 2, a ^ 0) as 

e ^ a-2/3$i((7-^/3(E - 2)), p ^ (7-2/3^2 K^/'(^ - 2)), (3.2) 

where the scaling functions $i and $2 are known analytically. 

Roughly speaking, the eigenstates whose energy lies near the band edges have a lo- 
calisation length of order a"^/*^. These states are therefore much more localised than 
typical eigenstates within the band, whose localisation length is of order ^0 ~ cr~^- Only a 
small fraction of the whole spectrum, of order u^/^ or ^q^^^, consists of these anomalously 
localised eigenstates. 

3.2 Heuristic analysis of the localised regime 

We now turn to a heuristic investigation of the moments of the position and of the 
probability density, in the presence of a weak diagonal disorder, and in the localised regime 

(t>eo>i). 

Consider for definiteness the Anderson model on a very long chain made of N :$> 1 
sites, for a given realisation of the random potentials {vn}- Let E"' be the energy eigen- 
values, labelled in some way by an integer a, and i/j^ be the corresponding eigenvectors. 
We have 

E "^n^n = 5"'^ E = Sm,n. (3.3) 

n a 
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We define the centre-of-mass co-ordinate n^^ and the locahsation length ^'^ of every eigen- 
state as 

<m = E^«)'' r = (E(^ (3.4) 

n n 

The initial condition (1.2) can be expanded as 'i/'n(0) = Sn,o = Z^aV'nV'o- ^ 
consequence, the wavefunction reads at all times t > 

^n(t) = $^e-'^"*CV'o"- (3.5) 

a 

3. 2. 1 Moments of the position 

Let us first take the example of the mean squared position, for which eq. (3.5) yields 

M2{t) = Y.e-^^^"-^'H^4Y.^^rni^'n. (3.6) 

a,/3 n 

Our heuristic analysis of this expression will be based on the following two hypotheses. 

(A) Interference terms between different quantum states can be neglected for large enough 
times. Eq. (3.6) thus becomes in the localised regime 

M2(oo)^5](C)'E"'(0'- (3.7) 

a n 

(B) Scaling properties of eigenstates can be modelled by considering that the probability 
density (V'n)^ roughly uniform over the range \n — n'^^\ < This simple scaling hy- 
pothesis has been shown by analytical means to hold in a variety of models, including ran- 
dom band matrices [7, 8], and the continuum Schrodinger equation in one dimension [13]. 
It amounts to stating that single eigenstates of the one-dimensional Anderson model do 
not exhibit multifractality, in contrast with earlier claims based on numerical evidence [22] . 

Consider first an eigenstate localised near the origin (|n"jjj| ^ ^"). For such an 
eigenstate, we have Xln^^(V'n)^ ~ (^")^5 while the prcfactor (i/Jq)'^ scales as 1/^". Now, 
in a small energy interval AE around some energy E, there are altogether Np(E)AE 
eigenstates, among which only a finite number, of order ^{E)p{E)AE, have \n^^\ ~ 

All these eigenstates bring comparable contributions, of order {^{E)) , to the sum in 
eq. (3.7). It is worth noticing that the factor $,{E) in the number of relevant eigenstates 
just compensates the prefactor (V'o )^ ~ We are thus left with the estimate 

M2(oo)~((r)'>„, (3.8) 

where the angular brackets denote an average over the whole spectrum of eigenstates a. 
More explicitly, 

M2(oo)- / {aE))%{E)dE^Cl (3.9) 
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We thus obtain the physicaUy intuitive result that the mean squared position saturates to 
a value of order for t ^ ^q. Generalising the above argument, we get the asymptotic 
result 

M,(oo)^Ag^o^ (3.10) 

for all the moments of the position, deep in the localised regime. 

This prediction can actually be made quantitative, using results from the Russian 
literature [14-17]. The long-time density correlation function has been calculated in these 
references, for the continuum Schrodinger equation with a weak white-noise potential. 
If the initial wave packet is peaked in energy around some mean E, one has Mq{oo) ^ 
Pq(2$,{E)Y , with the notation of ref. [16], where the amplitudes Pq have been calculated 
analytically, for integer values of q. In the present situation, by averaging this prediction 
over the whole spectrum of energies, we recover eq. (3.10), with 

Aq=pq J^^ {2 sin' py. (3.11) 

We have in particular 

3C(3) , 7(l80C(5)+7r'^) 
A2 = = 0.901543, ^4 = -^^ = 15.5343, etc., (3.12) 

where ( denotes Riemann's zeta function. 

3.2.2 Moments of the probability density (participation numbers) 

Let us now turn to the more interesting case of the moments Sq{t) of the probability 
density, with q = 2,3, . . . being an integer. Eq. (3.5) and hypothesis (A) yield 

Sq{oc)^ J2 W)'---(V'o')'E(^nO'---(0'- (3.13) 

cei,...,aq n 

Let us assume that the eigenstates ai, . . . , are ordered according to increasing local- 
isation lengths: < < ... < ^"1, and employ again hypothesis (B). If all the q 
eigenstates have their centre-of-mass co-ordinates close enough to the origin (|n"^| <C ^""^ 
foTk=l,...,q), the product (V'^i)^ . . . (V'n')^ is nonzero where the q eigenfunctions have a 
good common overlap. This occurs in the intersection of all their ranges, i.e., for |n| < 
where the product is of order 1/ (^"^^"^ . . .^"''), while the sum over n brings a factor of 
^"^i . Here again, the number of relevant eigenstates in some energy range AE cancels out 
with the prefactors (•00 0^ • • • (V'o')^' whence the estimate 



\ S • • • S / ai,...,aq 

The averaging of this expression over the eigenstates cti, . . . , ctg is more subtle than in 
the case of eq. (3.8), since negative powers of the localisation lengths are involved. Hence 
the anomalously localised eigenstates near the band edges can, and indeed will, play a role. 
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• If all the eigenstates cti, . . . , ct^ belong to the bulk of the spectrum, their localisation 
lengths scale as and we obtain the normal prediction Sq{oo) ~ ^o^'^~^^ ■ 

• If, on the contrary, some of the eigenstates, namely ai, . . . ,am, with m > 1, belong 
to the band edges, while Om+i, • • • , ctq belong to the bulk of the spectrum, the quantity 
to be averaged now scales as ^^^^ ^'^ while the total fraction of such g-uples of 
eigenstates is of order The optimal choice of the number m of anomalous eigenstates 
is m — hence the anomalous estimate 5*^(00) ~ '■^'^ 

The exponents of the above two estimates coincide with those obtained in section 2.3 
in the absence of disorder, with the localisation length scale ^0 replacing time, the bulk 
of the spectrum replacing the allowed region, and the band edges replacing the ballistic 
fronts. Generalising the above reasoning to non-integer values of the index q, we predict 
the power-law behaviour 

S,{^)^B,^-^^''^ (Co»l) (3.15) 

for the participation numbers in the localised regime, with the exponent r(g) being given 
in eq. (2.20) and in Table 1. Finally, in analogy with the result (2.29) in the absence of 
disorder, a logarithmic correction of the form 

S,ioo) ^ hhk±J^ » 1) (3.16) 

is expected in the marginal case (g = 2) of the usual participation number. The values of 
the prefactors Bg, and A and /i cannot be predicted by this heuristic analysis. 

3.3 Scaling laws in the crossover regime 

So far, we have obtained two kinds of predictions concerning the moments Mq{t) of the 
position of the particle and Sq{t) of the probability density in the weak-disorder regime. 
On the one hand, the analytical results (2.17), (2.19) obtained in the absence of disorder 
are expected to hold more generally in the ballistic regime, i.e., for 1 ^ t ^ ^o- On 
the other hand, a heuristic scaling analysis led us to the predictions (3.10), (3.15) in the 
localised regime, i.e., for 1 <^ ^0 <^ t. 

It turns out that the exponents involved in these results always match between the 
ballistic and the localised regime, both in the normal and in the anomalous case. We are 
thus led to conjecture that the crossover between these two limiting situations is described 
by universal scaling functions of the single variable 

x=^, (3.17) 

throughout the scaling region where t and ^0 cire simultaneously large, with the ballistic 
regime corresponding to x <^ 1, and the localised regime to x 1. It is worthwhile to 
recall that the absence of an intermediate diffusive regime between the ballistic one and 
the localised one is a peculiarity of the one-dimensional geometry. 
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3. 3. 1 Moments of the position 

We thus propose the following one- variable scaling law for the moments of the position 
of the particle: 

Mq{t)^aqt''Fq{x), (3.18) 

where the amplitudes are as in eq. (2.18). The result (2.17) in the absence of disorder 
is recovered as -Fq(O) = 1, while the estimate (3.10) in the localised regime implies the 
power-law fall-off 

i^g(a;) ~ — a;-'? (a; > 1). (3.19) 

(Xq 

The scaling law (3.18) can be shown to hold as x ^ 1, by means of a direct perturbative 
expansion in the random potentials [23] . This approach, following the lines of refs. [24, 3] , 
yields 

Fq{x) = \-a^^'^x + --- (a;<l), (3.20) 
at least when the index q = 2k is an even integer, and in particular 

= — = 3.39531, = = 2.86715, etc. (3.21) 

37r 1357r 



The analysis of the continuum Schrodinger equation [17] yields a behaviour similar 
to cq. (3.20) at small ,t, as well as a singular correction of relative order (ln,T)/,T to the 
leading power law (3.19) at large x. This logarithmic correction has been included in the 
analysis of numerical data on large random band matrices [9]. 

3.3.2 Moments of the probability density (participation numbers) 

Similarly, we postulate the following one-variable scaling laws for the moments of the 
probability density 

Sg{t)^bqt-^^'i^Gq{x) (g^2), (3.22) 

where the amplitudes bq are as in eq. (2.21) in the normal regime {q < 2), and as in 
eq. (2.22) in the anomalous regime {q > 2). The result (2.19) is recovered as Gq{0) = 1, 
while the estimate (3.15) implies the power-law behaviour 

Gqix) ^ ^ x^^'i^ (x>l). (3.23) 

bq 



In the marginal case of the usual participation number S2{t), we expect a logarithmic 
correction of the form 

The result (2.29) in the absence of disorder is recovered as 

0(0) = ^ = 0.050660, x(0) = ^ ^^27^ = 0.239933, (3.25) 
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while the estimate (3.16) in the locahsed regime implies 



(f){x) Xx, x(a;) fti (// — Alna;)a; (x » 1). (3.26) 

3.4 Numerical results 

In order to confirm the scaling predictions of section 3.3, we have performed direct 
numerical simulations of the dynamics of the tight-binding Anderson model. Introducing 
a finite time step s, we have discretised eq. (1.1) into the difference equation 

lljn{t + e)= Mt - - 2i£(V'n+l(^) + i^n-lit) + ^^nV'n(^))- (3.27) 

In the absence of disorder, the dispersion relation of the corresponding stationary equation 
between energy E and momentum p now reads 

sin(£^;) = 2£cosp. (3.28) 

For a fixed momentum p, eq. (3.28) has a first solution E^^^ = 2 cosp + (4£^/3) cos^ p + ■ ■ ■ 
that smoothly converges toward the continuous-time expression (2.2), and a second solution 
£■(2) = Tr/g — £■(!) rnodulo 2n/e, corresponding to fast oscillations over a time scale e. 

We have taken the initial condition (1.2) at time t = 0, while we have chosen at time 
t = e the Taylor expansion of the solution of eq. (1.1), to second order in e included. The 
only non-zero components in the full initial condition read 

MO) = 1, Me) = 1 - me - (1 + vl/2)e^ 
^±i(£) = -ie - {vo + v±i)e^/2, V'±2(£) = -e'/S- 

This prescription reduces to the level of the amplitude of the fast oscillations related to 
the E^"^"^ branch of the dispersion relation (3.28). 

We want to emphasise that the discrete-time difference equation (3.27) has the very 
same physical contents as the continuous-time differential equation (1.1). As an illustration 

of this point, the mean squared position M2^\t) still obeys the ballistic law (2.14), albeit 
with a prefactor given by 



/ ^£^f!^ = l^il^!£!)^ = 2 + 2.V... (3.30) 

B 27r 1 — 4£^ cos^ p 

The value a2 = 2, characteristic of the continuous-time equation [see eqs. (2.14), (2.16)], 
is thus recovered, with a small correction in . 

We have performed numerical simulations of the difference equation (3.27), with a 
time step e — 0.05, so that the oscillations and other discretisation effects are negligible. 
The random site potentials have been drawn from a uniform distribution over the interval 
-Wjl <Vn< W/2, so that = W^/12 and 

Co - (3.31) 



14 



The data to be presented below correspond to the following ten values of the localisation 
length scale: 

^0 = 25, 35, 50, 70, 100, 140, 200, 280, 400, 560. (3.32) 

For each value of the strength of disorder W is taken from the pcrturbative rela- 
tion (3.31). The measured quantities are averaged over 1,000 independent realisations 
of the random potentials. For each realisation, eq. (3.27) is integrated up to a time 
^max = 5^0: Order to enter into the localised regime, and for a range of space \n\ < rimax, 
with nmax = 2.5 tmax = 12.5 ^Q, in order to fully encompass the ballistic fronts. 

3.4-1 Moments of the position 

We have first checked the validity of the scaling law (3.18) for the moments Mq{t) of 
the position of the particle, as well as the analytical predictions concerning the associated 
scaling functions, on the examples of q — 2 and q — 4. Figures 2 and 3 respectively show 
log-log plots of numerical data for M2{t) and M4{t), divided by their expressions (2.14) in 
the absence of disorder, against the scaling variable x — t/^Q. The data, corresponding to 
the six largest values of eq. (3.32) for ^q? collapse in a nice way. This demonstrates the 
existence of the one- variable scaling functions F2{x) and ^4(2;). 

The dashed lines show the linear correction (3.20), (3.21) of the scaling functions at 
small as well as their asymptotic behaviour (3.19) at large x, namely ^2(0;) ~ 0.4508/x^ 
and F4^{x) ^ 2.589/a;^. The full lines show the one-parameter phenomenological fits 
F2{x) = (l + 3.395a; + 0.154xln(x + 1) + 2.218a;2)-i and ^4(0;) = (1 + 1.4336a; + 
0.622 a;ln(a; + 1) + 0.6215 a;^)-^. These expressions incorporate the small-x and large-x 
behaviour recalled just above. The fitted parameters are the amplitudes of the x ln{x + 1) 
terms, refiecting the structure in (In a;)/a; of the leading correction term at large x discussed 
below eq. (3.21). The good quality of the fits shows the quantitative agreement between 
analytical predictions and numerical data. We have also checked that the amplitudes of 
the behaviour (3.19), (3.20) of the scaling functions at small and large x are recovered 
within better than 10 percent if they are left as free fitting parameters, instead of being 
imposed as constraints. 

3.4-2 Moments of the probability density (participation numbers) 

We now turn to the moments Sq{t) of the probability density, for which fewer analytical 
predictions are available. We have first investigated the anomalous scaling laws of these 
moments in the localised regime, on the example of q — 3. In order to check the power 
law (3.15), as well as the leading correction to it, which can be expected to be of relative 

order Co"^^^? "^^ have plotted in Figure 4 the product ^053(00) against ^l^^- The data 
points correspond to all the values of eq. (3.32) for ^o- For each the data for Ss{t) 
in the range 1 < t/^o < 5 have been extrapolated, in order to get a reliable estimate for 
5*3(00). The error bars on the numbers obtained in this way are comparable to the symbol 
size. The plotted data nicely follow the least-square fit y — 2.46 a; — 2.21, confirming thus 
both the power-law behaviour (3.15) and the anticipated nature of the correction term, 
and yielding the estimate 

B3 ^ 2.5. (3.33) 
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In order to check the logarithmic behaviour (3.16) of the participation number in the 

locahsed regime, we have plotted in Figure 5 the product S2{oo) against In^o- The data 
points have been obtained as those of Figure 4. They nicely follow the least-square fit 
y = 0.253 a; + 0.642, confirming thus the behaviour (3.16), and yielding the estimates 

A !^ 0.25, fi ^ 0.64. (3.34) 

We have then determined the full one- variable scaling functions, defined in eq. (3.22), 
for the participation numbers throughout the crossover from the ballistic to the localised 
regime. We have again considered the example of = 3. Figure 6 shows a log-log plot 
of S2{t), divided by its behaviour (2.19), (2.23) in the absence of disorder, against the 
scaling variable x = t/^Q. The data, corresponding to the six largest values of eq. (3.32) 
for again collapse in a nice way, demonstrating the existence of the scaling function 
G3{x). The asymptotic largc-x behaviour Gs{x) ~ 33.6 x^/"^, shown as a dashed line, is 
accurately obeyed for values of the scaling variable as small as a; ~ 0.5. The full line shows 
the one-parameter phenomenological fit G3(x) = (1 — 0.56a; -|- 67.8 a;^)^/^, with the fitted 
parameter being the amplitude of the middle term. 

We end up with an investigation of the logarithmic behaviour (3.24) of the participa- 
tion number (t), for generic values of the scaling variable x. To do so, for any fixed value 
of the ratio x = t/^o, we have performed a least-square fit of all the available data for tS2{t) 
against Int, the slope and the intercept of those fits respectively yielding estimates for (p{x) 
and x(^)- Figure 7 shows a log-log plot of the amplitude (t){x) thus obtained, against the 
scaling variable x. Dashed lines show the value of (j){0) given in eq. (3.25), and the asymp- 
totic behaviour (f){x) ~ 0.253 a; [see eq. (3.26)]. The full line shows the one-parameter 
phenomenological fit (f){x) = (0.00257 0.00282x + 0.0639x^)^/2, ^j^j^ ^^ic fitted parame- 
ter being again the amplitude of the middle term. Similar, albeit less accurate, numerical 
results have been obtained for the function x(^)- 

4 DISCUSSION 

The most salient outcomes of this work are summarised in Table 1. In the absence 
of disorder, and for an initially localised wave packet, the probability density of a tight- 
binding particle exhibits anomalous scaling and bifractality. Indeed its moments of order 
q, namely the dynamical participation numbers Sq{t), scale in time with a non-trivial 
exponent T{q) for q > 2. 

To put it more boldly, a free quantum-mechanical particle is bifractal. This striking 
feature is entirely due to the presence of ballistic fronts. For the tight-binding model 
considered in this paper, these fronts correspond to the transition region in the theory of 
Bessel functions. Both the existence of ballistic fronts and their width scaling as t^^^ are 
actually general characteristics of difference equations, with a bounded energy dispersion 
curve. To the best of our knowledge, this bifractality phenomenon has been overlooked so 
far. Plots similar to our Figure 1, showing the squared Bessel function with its ballistic 
fronts, have been displayed and described e.g. in ref. [25] , but without the authors noticing 
the relevance of the fronts. The moments Mq{t) of the position of the particle are not 
affected at all by the presence of ballistic fronts. Ref. [5] contains a rigorous and general 
discussion on the spreading of a wave packet. 
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For the Anderson model with a weak diagonal disorder, the bifractal phenomenon per- 
sists throughout the different regimes of the localisation phenomenon. Indeed anomalous 
scaling of the participation numbers, with the same exponent T{q), holds in the ballistic 
regime {t <^ ^o), in the localised regime {t 3> ^o), and in the crossover between them 
{t ~ Co), where these quantities obey scaling laws involving the single variable x = t/^Q. 
These scaling predictions have been confirmed quantitatively by accurate numerical simu- 
lations. 

The bifractal exponent T{q) characterises both the time decay of participation numbers 
in the ballistic case, and their dependence on in the localised regime. In the ballistic 
regime {t -C ^o): this phenomenon is related to the existence of a finite upper band edge, 
and to the uniform scaling of single eigenstates. This absence of multifractality is in 
tiun related to the direct crossover between a ballistic and a localised regime in the one- 
dimensional geometry, without an intermediate diffusive phase. These features have been 
established by means of a variety of analytical techniques [7, 8, 13-18]. In the localised 
regime {t S> ^o): bifractality is intimately related to the Derrida- Gardner anomalous scaling 
of the localisation length and of the density of states near band edges. 

Two very different physical mechanisms seem therefore to be responsible for the bifrac- 
tal behaviour observed in these two different regimes. Let us underline that both phenom- 
ena have nevertheless one important characteristic in common. They take place near the 
band edges of the energy spectrum (2.2) in the absence of disorder, where the two plane 
waves exp(±inp) become identical. This degeneracy of the plane-wave basis explains, at 
least qualitatively, both the fiatness of the dispersion curve around the upper band edge 
(hence the existence of sharp ballistic fronts) and the high sensitivity of eigenstates to a 
perturbation (hence the Derrida- Gardner anomalous scaling). 

The usual participation number 5*2 (t) corresponds to the borderline case {q = 2) 
between the normal and the anomalous regime. This quantity exhibits a logarithmic 
correction, either in time [see eq. (2.29)] or in [see eq. (3.16)]. A similar logarithmic 
behaviour has been described [26] in the case of the mean return probability 

C{t) = l ['poiu)du. (4.1) 
t Jo 

In the absence of disorder, this quantity has a logarithmic correction to the naive scaling 
C^^\t) ~ 1/t, which can be easily derived by the Mellin approach, yielding 

= i £ (J„(2„))^d„ « + (4.2) 

a result very similar to eq. (2.29). This behaviour was shown in ref. [26] to explain the 
occurrence at several places in the literature of a fake non-trivial decay exponent 6 ~ 0.84 
for the mean return probability, even in the absence of disorder, as well as a variety of 
erroneous conclusions drawn from there. 

Beyond the amusing phenomenon of bifractality demonstrated in the one-dimensional 
Anderson model, the present work underlines that the moments of the position of a 
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quantum-mechanical particle and the moments of the associated probability density (par- 
ticipation numbers) can exhibit different scaling laws, with unrelated dynamical exponents. 
This outcome corroborates the recent general discussion [5] on the complexity of the scaling 
laws governing the dynamics of quantum systems. 
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CAPTIONS OF TABLE AND FIGURES 



Table 1: Summary of various characteristic features of the bifractahty phenomenon. 

Figure 1: Plot of the probability density Pn^\t) = (j^(2t))^ in the absence of disorder, 
against site number n, for a time t — 100. Arrows indicate the semi-classical ballistic fronts 
at n = ±200. 

Figure 2: Log-log plot of the ratio of mean squared position M2{t) to its value M2^\t) 
in the absence of disorder, against the scaling variable x = t/^o- Symbols: numerical data 
for various values of ^o- Phenomenological fit (full line) and asymptotic behaviour (dashed 
lines) are given in the text. 

Figure 3: Same as Figure 2, for the fourth moment M4{t) of the position. 

Figure 4: Plot of the product -53(00) against Co^^- Symbols: numerical data for various 
values of ^o- The least-square fit (full line) is given in the text. 

Figure 5: Same as Figure 4, for the product ^0'S'2(oo) against In^o- 

Figure 6: Log-log plot of the ratio of the participation number 83(1) to its value S^^\t) 
in the absence of disorder, against the scaling variable x = t/^Q. Symbols and lines: as in 
Figure 2. 

Figure 7: Log-log plot of the amplitude (^{x) of the participation number S2{t), defined 
in eq. (3.24), against the scaling variable x = t/^Q. Symbols: numerical data. Lines: as in 
Figure 2. 
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Table 1 



component of bifractality 


normal 


anomalous 


relevant region of wavefunction 
in ballistic regime (t <^ ^o) 


bulk 

(allowed region) 


fronts 
(transition region) 

p (+\ ^ +-2/3 


relevant eigenstates 
in localised regime (t ^ Co) 


bulk of spectrum 
^ ~ ^0 ~ 


band edges 
2/3 


range of index q 


q<2 


g > 2 


exponent r{q) 


q-1 


2g- 1 
3 


Renyi dimension Dq 


1 


2q-\ 
3(g-l) 


multifractal spectrum 




q; = -, t = - 
3' 3 
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In (M^(t)/Mf (t)) 




In (S3(t)/sf (t)) 




In 0(x) 




